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Gravitational waves from the axial oscillation of neutron star 
considering non-Newtonian gravity 
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Abstract The eigen-frequencies of the axial w-modes of neutron star described by a super-soft equation of state 
(EOS) are investigated, by considering the non-Newtonian gravity. The results show that at the same stellar mass, the 
frequencies of w; and wp for our model are lower than that of the typical EOSs (such as APR); and the frequencies 
increase with the stellar masses, which is contrary to that of the typical EOSs. These characters may provide a probe to 
testify the super soft symmetry energy and the non-Newtonian gravity in the future. Moreover, our model also has the 
universal behavior of the mass-scaled eigen-frequencies as a function of the compactness. 
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Introduction 


certain 


supra-saturation densities" "^", The latter kind of 


Neutron stars (NSs) provide a perfect model for 
studying the gravitational waves and the equation of 
state (EOS) of the dense matters? 
important topics in modern physics. For both the 


l which are two 


topics, it is believed that there is a very useful type of 
oscillation mode, the w-mode, being an important 
mechanism for producing gravitational waves!) 
which will contain abundant information for the EOS 


[7-10] The w-mode only 


of the super dense matters 
exists in general relativity, which is associated with the 
space-time and for which the motion of the fluid is 
13A, related to the 


space-time curvature and exists for all relativistic stars, 


negligible The w-mode is 
including black holes. The standard axial w-mode is 
categorized as wy. 

At present time, however, the EOS of matter 
under extreme conditions is still rather uncertain. One 
of the main sources of uncertainties in the EOS of 


neutron-rich matter is the poorly known density 


dependence of the nuclear symmetry energy, Esym(p)!”. 


Many of the theories predict that the symmetry energy 


[12-16 


increases continuously at all densities ! But many 


other models predict that the symmetry energy 


symmetry energy functions is generally regarded as 
being soft. Among the soft models, the original 
Gogny-Hartree-Fock (GHF) model! and group Tr?! 
are super soft, in which the pressure drops quickly to 
zero around three times of the saturation density po. 
The super-soft models either cannot keep the NSs 
stable or predict maximum NS masses significantly 
below1.4M; (M. is the solar mass). On the other 
hand, the recent experimental results support the 
super-soft prediction by analyzing the FOPI/GSI 
experimental data on the z /z' ratio in relativistic 


B2] within a transport model by 


heavy-ion collisions 
using the MDI (Momentum-Dependent-Interaction) 
EOSP?l. In order to let the super-soft EOS support the 
mass observation of neutron star, the non-Newtonian 
gravity was introduced". In this work, the w-mode is 
investigated by using a super-soft EOS, where the 
non-Newtonian gravity is considered. 

In Section 2, we describe the EOS and the non- 
Newtonian gravity. In Section 3, we review the 
formalism used to calculate w-modes and present the 


results before concluding in Section 4. 
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2 Equation of state and non-Newtonian 
gravity 


The MDI is a phenomenological effective interaction 
based modified — finite-range Gogny 
interaction . By using the MDI EOS, the 
circumstantial evidence for a super-soft symmetry 


on a 
[26,33,35-37] 


energy was found recently from analyzing the 
FOPI/GSI experimental data"?! on the r/r ratio in 
relativistic heavy-ion collisions"), For a neutron star 
consisting of neutrons, protons, and electrons (npe), 
the internal pressure of neutron stars is given by 
Eq.(1)83437) with the chemical equilibrium conditions 
of Me=Un-Mp=4 0E,(p) and charge neutrality requirement 
of p=(1-0)p/2, 


dE, dE,(p) 
Paed =| o 4 2 


lilii () 


*z00-3)pE, (o). 


The energy nucleon for symmetric nuclear 


matter can be well approximated as P^ 38] 


8x B " 
Ep) = gy Py *45-C216554— (7) 


5mh! p Po o+l p 
1 A y a2 2 2 2 

+ C, 4 C, A^x 6p; -A 

dno" 1 WG? [p;C Py ) (2) 


2p 
—8Ap! arctan ——- 
D, A 


Ts oan, 4p +A? 
uU +1247 pm, 


and the symmetry energy can be expressed ag 9l 


4Bx Bx z 
] (P 


8 s, p 
E = + 24.594 
(P) Omit p ^! 4p, otl o+l p, 


4p; * A (3) 
A? ] 


C, ,Am i2 4 2:2 
—) A‘[4p,-A In 
Spp E [4p, - ^ p, 


LG 
95,p 


4r Ap; +h? 
Ga) Tp; — p; Ap; * A) 5] 


5 is the Fermi momentum for 


where p, -h(3z^ p[2) 
symmetric nuclear matter, and the values of the 
parameters are o=4/3, B=106.35 MeV, Cy—11.70 
MeV, C,=-103.40 MeV and A=p)=p,(p), 
where o, is the saturation density. The value of Eyym(p) 
is 31 MeVP* 9 The parameter x was introduced to 
vary the density dependence of the symmetry energy 
without changing any property of symmetric nuclear 
matter. It has been shown that x=1 can reproduce the 
FOPI/GSI pion production dataP?! within the transport 


model analysis, and corresponds a  super-soft 
symmetry energy. The corresponding EOS is denoted 
as MDIxI. 

In the traditional framework, MDIx1 EOS 
cannot support the observed masses of pulsars because 
of its super-soft symmetry energy. To support the 
Observations of neutron stars, the non-Newtonian 
gravity was introduced". The non-Newtonian gravity 
can be described by adding a Yukawa term to the 
conventional gravitational potential between two 


objects 1. that 1s, 


V(r) = Bu * ae ^), (4) 
r 


where G is the universal gravitational constant, a is a 
dimensionless strength parameter, and A is a length 
scale. In the scalar/vector boson exchange picture 
a=+g"/(4nGmy_) and A=1/u, where g, u and m; are the 
boson-baryon coupling constant, the boson and baryon 
Within the 


non-Newtonian gravitational 


mass, respectively. mean-field 


approximation, the 
effects can be described by the equation of state of 


dense matters, but Einstein equation remains 


unchanged. The extra energy density due to the 


non-Newtonian gravity can be calculated by 


where V is the normalization volume, p is the baryon 
number density, and 7=| x,—x, |. The corresponding 


addition to the pressure is then 
Pos = Eug -gplQu) (6) 


where a constant boson mass independent of the 
density is assumed" ^^! Including the non-Newtonian 
gravity the EOS becomes P=Prpet+Pup, where Prpe is 
the conventional pressure for neutron star mainly 
consisting of neutrons, protons and electrons. Shown 
EOS of MDIxI, 
non-Newtonian gravitational parameter gh? is taken 
as 75 GeV? or 100 GeV? (denoted as MDIx1-75, 
MDIx1-100). In order to show the effect of the 
non-Newtonian gravity, the typical EOS, APR is 
used for comparison. The mass-radius relation is 
displayed in Fig.2. It is easy to see that for the EOS 
MDIkx1, the radius of the neutron star is far larger than 
that of the typical EOS APR. 


in Fig.l is the where the 
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Fig.1 The equation of states versus the reduced baryon 
number density p/p». 
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Fig.2 The mass-radius relations. 


3 The gravitational waves from the axial 
w-modes of neutron stars with super soft 
EOS and considering non-Newtonian gravity 


First, we briefly introduce the formalism for the axial 
w-modes. Here we use geometrical unit (G=c=1). 
According to Chandrasekhar & Ferrari"l, the axial 
perturbation equations for a static neutron star can be 
simplified by introducing a function z(r), constructed 
from the radial part of the perturbed axial metric 
components. It satisfies the following differential 
equation 

dz 5 

— -* [o -V(r)]z - 0, (7) 

dr? 
where o (=@+iq@;) is the complex eigen-frequency of 
the axial w-mode, r, is the tortoise coordinate"). 
Inside the star, the potential function V is defined by 


- [4 Dr Azr^ (p p) - 6m] (8) 
F 


V = 


where / is the spherical harmonics index, p is the 
density, p is the pressure, m is the mass inside radius r, 
and e" and e^ are the metric functions given by the line 
element for a static neutron star. Outside the neutron 
star, Eq.(8) reduces to 

r-2M 


4 
r 


V= [I -- Dr - 6M], (9) 


where M is the total gravitational mass of neutron star. 
The solutions to this problem are subject to a set of 
(BC) 


and Ferrari The 


boundary conditions constructed by 
Chandrasekhar 
calculations are performed by applying the continued 
fraction method" ^?! together with the EOSs described 


in Section 2. 


numerical 


v/ kHz 


10 12 14 16 
M | Mo 


18 2.0 


Fig.3 Frequency (a) and damping time (b) of w;-mode as a 
function of the neutron star mass M. 

In what follows we present our numerical 
results for the axial w-modes of neutron stars 
described by our models. Shown in Fig.3 are the 
frequencies (a) and the damping times (b) of the 
wrmode as a function of the stellar mass. In order to 
show the effect of the non-Newtonian gravity on the 
w-mode, we also plot the corresponding results of 
typical EOS, APR, which is often used as a comparing 
criterion in the neutron star studies. Fig.3(a) shows 
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that at the same stellar mass, the frequency of the 
wrmode in the neutron star described by MDIx1 with 
non-Newtonian gravity is smaller than that of the 
typical EOS APR. Moreover, for our model, the 
frequencies increase with the stellar mass, which is 
contrary to that of the typical EOS APR. The 
difference of the frequencies between our model and 
the typical EOSs may come from the neutron star 
globe structure, as shown in Fig.2, that is, the neutron 
stars described by our model have larger radii. And a 
larger non-Newtonian gravitational parameter glu 
corresponds with a small frequency. The damping time 
of wr-mode showing in Fig.3(b) indicates that there is 
no notable distinguishability among the three 
considered EOSs for the neutron star with canonical 
mass1.4M,. 


—«— MDIx1-75 
0.20 ; --&-- MDIx1-100 
--e-- APR 


LJ 
0.08 0.12 0.16 0.20 0.24 
MIR 


Fig.4 Frequency (a) and damping time (b) of w,-mode scaled 
by the stellar mass M as a function of compactness M/R. 
Previous work has shown that the eigen- 
frequency of w-mode scaled by the stellar mass 
exhibits a universal behavior independent of the EOS 
as a function of the compactness M/R, which allows an 
accurate determination of the eigen-frequency of the 
gravitational waves from the w-mode if the masses and 
radii of neutron stars are accurately observed"?! [n 
Fig.4, we show the frequency (a) and the damping 
time (b) of wr-mode scaled by the stellar mass M as a 
function of compactness M/R. Also, our results have 
the universal behavior, and the scaled frequencies of 
the two models described by MDIx1-75 and 
MDIx1-100 have a perfect universal behavior. 
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Fig.5 Frequency (a) and damping time (b) of wp-mode as a 
function of the neutron star mass M. 
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Fig.6 Frequency (a) and damping time (b) of wp-mode scaled 
by the stellar mass M as a function of compactness M/R. 

We also investigate the eigen-frequencies and 
their scaled properties of wp-mode. The frequencies (a) 
and the damping times (b) of the wp-mode as a 
function of the neutron star mass M are shown in Fig.5, 
and these quantities scaled by the stellar mass M as a 
function of compactness M/R are shown in Figs.6(a) 
and 6(b). It is shown in Fig.5(a) that the frequencies of 
Wp-mode described by the MDIxl1 are almost a 
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constant for a fixed non-Newtonian gravitational 
parameter, e.g., 12.6 kHz for g’/’=75 GeV? and 
10.9 kHz for g’/u’=100 GeV”. The frequencies of the 
APR model descend obviously as the stellar mass 
increases. Similarly, there also exists the scaled 
universal property for the wp-mode, as shown in Fig.6. 


4 Conclusion 


In this work, we calculated the eigen-frequencies of 
the axial w-modes of neutron star described by a 
super-soft equation of state (EOS), where the 
non-Newtonian gravity is considered. It is found that 
the at the same stellar mass, the frequencies of wj and 
Wp for our model are lower than that of the typical 
EOSs (such as APR). And the frequencies are 
increasing as the stellar mass increases, which is 
contrary to that of the typical EOSs. The characters of 
the frequencies of w;and wp can be used as a probe to 
testify the super soft symmetry energy and the 
non-Newtonian gravity in the future. Moreover, we 
also show that this model also has the previously 
behavior of the 


eigen-frequencies as a function of the compactness. 


found universal mass-scaled 
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